Astronomy & Astrophysics manuscript no. draft-astroph 


© ESO 2013 


February 13, 2013 





X 



The role of low-mass star clusters in massive star formation. The 

Orion Case. 

V. M. Rivilla 1 , J. Martm-Pintado 1 , I. Jimenez-Serra 2 , and A. Rodriguez-Franco 1 

1 Centro de Astrobiologfa (CSIC-INTA), Ctra. de Torrejon-Ajalvir, km. 4, E-28850 Torrejon de Ardoz, Madrid, Spain 
e-mail: rivilla@cab.inta-csic.es 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., 02138, Cambridge, MA, USA 



m 
o 

(N 



< 

o 

6 



> 

m 

(N 

o 



Received 16/06/2011; accepted 28/01/2013 



ABSTRACT 



Context. Different theories have been proposed to explain the formation of massive stars: two are based on accretion 
processes (monolithic core accretion and competitive accretion), and another on coalescence of low- and intermediate- 
mass stars. To distinguish between these theories, it is crucial to establish the distribution, the extinction, and the 
density of low-mass stars in massive star-forming regions. 

Aims. Our aim is to analyze deep X-ray observations of the Orion massive star-forming region using the Chandra Orion 
Ultradeep Project (COUP) catalog, to reveal the distribution of the population and clustering of PMS low-mass stars, 
and to study their possible role in massive star formation. 

Methods. We studied the distribution of PMS low-mass stars with X-ray emission in Orion as a function of extinction 
with two different methods: a spatial gridding and a close-neighbors method, with cells of ~ 0.03 x 0.03 pc 2 , the typical 
size of protostellar cores. We derived density maps of the stellar population and calculated cluster stellar densities. 
Results. We found that PMS low-mass stars cluster toward the three massive star-forming regions: the Trapezium 
Cluster (TC), the Orion Hot Core (OHC), and the OMC1-S region. We derived PMS low-mass stellar densities of 10 5 
stars pc -3 in the TC and OMC1-S, and of 10 6 stars pc -3 in the OHC. The close association between the low-mass 
star clusters with massive star cradles supports the role of these clusters in the formation of massive stars. The X-ray 
observations show for the first time in the TC that low-mass stars with intermediate extinction are clustered toward 
the position of the most massive star 6 Ori C, which is surrounded by a ring of non-extincted PMS low-mass stars. 
This 'envelope-core' structure is also supported by infrared and optical observations. Our analysis suggests that at least 
two basic ingredients are needed in massive star formation: the presence of dense gas and a cluster of low-mass stars. 
The scenario that better explains our findings assumes high fragmentation in the parental core, accretion at subcore 
scales that forms a low-mass stellar cluster, and subsequent competitive accretion. Finally, although coalescence does 
not seem a common mechanism for building up massive stars, we show that a single stellar merger may have occurred 
in the evolution of the OHC cluster, favored by the presence of disks, binaries, and gas accretion. 
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1. Introduction 

Massive stars (M>8M Q ) are the main source of mechanical 
energy injection into the interstellar medium of galaxies. 
Although they are important for the evolution of galaxies, 
their formation process is not clearly understood. This is 
due to the large distances to the massive star-forming re- 
gions, their short evolution time-scales, the complexity of 
their birth regions and the high extinction. Several theories 
have been proposed to explain the formation of the most 
massive stars: i) monolithic gravitational collapse and core 
accretion (Yorke fc Son nhalter, l2"UfJ2l McKee & Tan, [217051 
Krumholz et al., 12009^ : ii) coalescence of low- mass stars 
in high stellar density clusters (Bally & Zinnecker, |2005 
Bonnell et al., [TM51 Stahler et al, IMJOj) ; and hi) com- 
petitive accretion of low- and intermediate-mass stars in 
the cluster potential well (Bonnell & Bate, I2006|) . While 
the classic theory of monolithic core accretion predicts the 
formation of mainly massive companions in small multi- 



ple systems, coalescence and competitive accretion models 
require a dense population of low-mass stars that form a 
cluster. Therefore, to distinguish between the different sce- 
narios of massive star formation, it is crucial to establish 
the distribution and density of low-mass stars in massive 
star cradles. 

Measuring the population of low-mass stars where the 
early phases of massive star formation take place has been 
extremely difficult because they lie behind high extinction. 
X-ray observations are particularly useful for studying high- 
mass star-forming regions because they can penetrate into 
the cloud despite the high extinctions found in these regions 
(>15 mag). In addition, they suffer less from foreground- 
and background galactic contamination than optical or in- 
frared (IR) studies, and allow a census of pre-main sequence 
(PMS) low-mass stars embedded in the molecular cloud. 

Orion is the nearest (d~414 pc) region of high-mass 
star formation, where three sites of massive star forma- 
tion have been identified: the Trapezium Cluster (TC), the 
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Orion Hot Core (OHC), and OMC1-S. The TC contains the 
well-known four main-sequence (MS) massive Trapezium 
stars and is located in the central region of the optically 
visible Orion Nebula Cluster (ONC). The ONC has a size 
of —3 pc (Hillenbrand & Hartmann, I1998P and an age of 
about 10 6 yr (Hillenbrand. HMT)) . The OHC, locat ed within 
the Orion KL nebula (Zuckerman et al., HIM Shuping 
et al., I2U04P is a very compact (size of ~ 0.03 pc), hot 
(T>200K), and dense molecular condensation (mean den- 
sity of n[H 2 ]=5x 10 7 cm -3 , Morris et aL. HTjgtj]), w ith a total 
luminosity of few 10 5 L (de Vicente et al. . 130113)) . The OHC 
is the prototype object for studying the formation of mas- 
sive stars. Three massive objects with masses around 10 M Q 
have been identified in this region: the Becklin-Neugebauer 
Object (BN), and sources / and n. 

OMC1-S ( ~10 4 L Q ; Mezger et al., [T39U1 Rodriguez- 
Franco et al., I1999ap . located ~ 90" south of the OHC, 
is the other young massive star-forming region embedded 
in the Orion molecular cloud. OMC1-S is a very active star 
forming region, harboring several outflows (Ziurys et al., 
1990; Rodriguez-Franco et al., 1999a, 1999b; Zapata et al., 
I3PP51 1501151) . 

The Chandra Orion Ultradeep Project (COUP; Getman 
et al.. 12005a) . based on a 10-day Chandra X-ray observa- 
tion centered on the ONC, provided a unique data base 
of X-ray PMS stars. The X-ray emission from these low- 
mass stars can mainly be attributed to plasma heated to 
high temperatures by violent reconnection events in mag- 
netic loops in the corona of PMS low-mass stars. Similar 
sources, variable in X-rays and in radio, have been detected 
in other low-mass star-forming regions with X-ray proper- 
ties consistent with coronal magnetic activity (Feigelson & 
Montmerle, [TM1 Forbrich et al.. |2H0T|) . 

In this paper, we extend previous works based on the 
COUP database to study the spatial distribution and esti- 
mate the densities of low-mass PMS stars. We find low-mass 
PMS stars clusters toward the three regions where massive 
star formation has recently occurred (TC), or it is currently 
taking place (OHC and OMC1-S). This fact and the high 
densities of PMS stars in these regions suggest that low- 
mass stars could play a role in the formation of massive 
stars. 

2. Data analysis 

We analyzed the distribution and stellar densities of PMS 
low-mass stars with X-ray emission in Orion as a function of 
extinction. We used the Chandra COUP database provided 
by Getman et al. (20f J5a|) . We applied two different meth- 
ods: i) spatial gridding, and ii) a close-neighbors method. 
This section is structured as follows: In 2.1. we present the 
two methods. Sections 2.2, 2.3, 2.4, and 2.5 are devoted 
to the effects introduced in our analysis by the cell size, 
the selected extinction ranges, the background/foreground 
contamination, and incompleteness. 

2.1. Methods 

2.1.1. Spatial gridding 

We integrated the spatial distribution of the X-ray stars 
into a grid of square cells (see section 2.2. about the cell 
size), counting the number of stars contained in each cell. 
This method is very adequate to obtain maps of the stel- 



lar density, because it samples the entire region. However, 
spatial gridding can underestimate the real values of the 
stellar density peaks (see section 2.1.2.). 



2.1.2. Close-neighbors method 

The clustering of low-mass PMS stars may be affected by 
the spatial gridding used to compute the density of X-rays 
sources across the Orion region. For instance, depending 
on the grid origin, a stellar cluster could be within a single 
cell, or could be divided into two or more adjacent cells, 
leading to different values of the peak stellar densities. To 
estimate this effect we also analyzed of the distribution of 
close neighbors centered on every COUP source, which is 
independent of the grid origin. This method allows us to 
calculate the real stellar density peaks, but it is unable to 
generate density maps because it undersamples the field. 



2.2. Cell size 

From IR observations, Lada et al. (|2004p found several 
dense stellar clusters in Orion, some of them with sizes of 
< 0.05 pc. We are interested in detecting of these compact 
clusters, and in particular those connected with massive 
star formation. Massive stars are born in dense and hot 
condensations of gas, with typical sizes in the range 0.01 
pc - 0.1 pc (Kurtz et al.. l2000p . with a mean size of 0.03 pc 
(Saito et al., 2008). This size also matches the size of the 
Orion Hot Core (Wright et al.. H996| ) . Therefore, we selected 
a cell size of 15" x 15" (0.03 pcx0.03 pc) to match the mean 
size of the molecular cores where massive stars form as well 
as the sizes of dense stellar clusters. To evaluate the effect 
introduced by the cell size, we also repeated our analysis 
using a larger cell of 24" x 24" (-0.05 pcx0.05 pc). We 
compare the results obtained with these two different cell 
sizes in section 3.2.1. 



2.3. Selected extinction ranges 

Because the extinction across the Orion central region 
changes by several orders of magnitude (from 10 21 to 10 24 
cm -2 ), we analyzed the X-ray population as a function of 
three X-ray extinction ranges. The selection of the extinc- 
tion bins are based on physical arguments. The first extinc- 
tion range, logN#< 22.0 cm~ 2 (Ay<5 mag), corresponds 
to the regions with low extinction that are mainly affected 
by the foreground gas in front of the optically visible ONC. 
With the bin logN#>22.5 cm -2 , corresponding to Ay>15 
mag, we expect to find sources that are heavily embedded 
in dense regions of hot gas and dust (hot cores), where it is 
believed that massive protostars are found. The condition 
logN#>22.5 cm" 2 has also been used by Wang et al. (|2007p 
to identify the heavily embedded population obscured by 
molecular material in the massive star-forming region NGC 
6357. The population with 22.0 cm" 2 <logN i/ <22.5 cm" 2 
(5<Ay<15 mag) corresponds to slightly extincted sources. 
In our analysis we only considered COUP sources with high 
enough X-ray counts to allow one to determine the value of 
logN# with an error of <0.5 dex. 
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Fig. 1. Spatial distribution of low-mass stars for three different extinction ranges: logN# <22.0 cm~ 2 (upper panels), 
22.0 cm _2 <logN#<22.5 cnr 2 (middle panels), and logN#>22.5 cm -2 (lower panels). Dashed contours represent dense 
gas traced by the integrated intensity emission of CN (N=l-0) (Rodriguez-Franco et al.. 11998]) . The first contour level 
corresponds to 7 K km s _1 , the interval between contours is 4 K km s _1 . The left panels show the position of the X-ray 
stars as blue dots. Four crosses, the open square, and the open triangle show the location of the four main sequence 
massive Trapezium stars, the OHC, and the OMC1-S. The Right panels report the stellar surface density for the same 
extinction ranges, derived by counting the number of COUP sources using a cell size of 15" x 15" (0.03 pc x 0.03 pc) 
(right color scale), superimposed on the CN emission. The color contours have been smoothed for display purposes. 
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Fig. 2. Zoom of the central region from Fig. 1. The color contours represent the stellar surface density (number of COUP 
sources using a cell box of 0.03 pcx0.03 pc, right color scale) for the three extinction ranges (logN# <22.0 cm -2 in the left 
panel; 22.0 cm -2 <logN#<22.5 cm -2 in the middle panel; and logN/f >22.5 cm -2 in the right panel). The color contours 
have been smoothed for display purposes. As in Fig. 1, the dashed contours represent the CN emission, and the four 
crosses, the open square, and the open triangle show the location of the four main-sequence massive Trapezium stars, 
the OHC, and the OMC1-S. 



2.4. Background/ foreground contamination 

One of the advantages of using X-rays to study massive 
star-forming regions is that they suffer less foreground- and 
background contamination than optical and infrared obser- 
vations. Galactic contamination has very little impact on 
X-ray studies, but cxtragalactic sources can still be con- 
fuse d with embedded PMS stars. However, as Getman et 
al. (|2005bp showed, the three regions studied in detail in 
this paper (i.e., the TC, OHC, and OMC1-S) are nearly 
free from contamination, so the detected sources reflect 
the actual cluster-member population. Despite this, we es- 
timated a non-zero background/foreground level to pro- 
vide a reasonable upper limit to the error of the stellar 
densities calculated in this work. A very conservative up- 
per limit for the contamination would be to consider a 
uniform distribution of all COUP sources (~1600) in the 
field observed by Chandra (17' x 17'). This gives an er- 
ror of ±0.35 stars cell -1 , considering 0.03 pcx0.03 pc cells. 
Getman et al. (2005b) claimed that only 217 COUP sources 
do not belong either to the ONC or the molecular cloud. 
Considering these sources, we obtain a lower limit for the 
background/foreground contamination of only 0.05 stars 
cell -1 . In our density analysis we consider a contamina- 
tion level between the two values of 0.2 stars cell -1 , which 
is equivalent to 3 x 10 2 star pc -2 in surface density and 2 
x 10 4 stars pc -3 in volume density. 

2.5. Incompleteness 

Although the COUP observations are quite deep and hence 
very sensitive, there is still a missing population of sources 
due to obscuration, especially in regions with high extinc- 
tion, like the OHC and OMC1-S. By comparing the X-ray 
luminosity functions of the incomplete OHC and OMC1- 
S populations with the more complete ONC population, 



Grosso et al. (2005) estimated that only about 48% of the 
low-mass sources have been detected in the OHC, and 63% 
in the OMC1-S region. The TC is not as strongly obscured, 
and therefore its population is considered to be nearly com- 
plete. 



3. Results 

3.1. Distribution of the low-mass stellar population 

We used spatial gridding to obtain maps of the stellar den- 
sity as a function of extinction. Fig. 1 shows the spatial 
distribution of X-ray low-mass stars for three different ex- 
tinction ranges (see section 2.1.) superimposed on the total 
column density of dense molecular gas traced by the emis- 
sion of CN (N=l-0) (dashed contours ; Rodriguez- Franco 
et al., I1998p . The left panels show the location of the X- 
ray stars as blue filled circles. Four crosses, an open square, 
and an open triangle indicate the location of the four TC 
stars, the OHC region, and the OMC1-S region. The right 
panels present the stellar surface density in terms of the 
number of COUP sources contained in cells of 15" x 15" 
(0.03 pcx0.03 pc). A zoom of the central part of the field 
is presented in Fig. 2. Consistent with Grosso et al. (2005), 
we found that PMS low-mass stars cluster toward the three 
regions related to massive stars: the TC, the OHC and the 
OMC1-S region. 

The less extincted population, with logN# < 22.0 cm -2 
(upper panels in Fig. 1, left panel in Fig. 2), shows a ring- 
like structure around the TC. The slightly extincted stel- 
lar population, with 22.0 cm -2 <logN#<22.5 cm -2 (mid- 
dle panels in Fig. 1 and 2), is clustered toward the posi- 
tion of the most massive star in the Trapezium, 9 1 Ori C. 
Remarkably, this population partially fills the hole found 
in the non-embedded population. The most extincted X- 
ray sources (logNjy>22.5 cm -2 ; lower panels in Fig 1., right 
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panel in Fig. 2), clearly associated with the molecular cloud 
traced by the CN emission, are clustered toward the OHC 
and the OMC1-S regions (open square and triangle, respec- 
tively, in Fig. 1 and 2), which host the ongoing massive star 
formation activity. 

The selected extinction ranges have major implications 
on the outcome of this type of studies. Previous studies 
by Feigelson et al. (|2HT75|) and Grosso et al. ([2T)tl5"[) only 
considered two extinction ranges with logN#>22.0 cm" 2 
and logN#<22.0 cm -2 . Although these authors already re- 
ported the clustering toward the TC, OHC, and OMC1-S, 
their classification criteria did not allow them to clearly sep- 
arate a non-extinctcd PMS population in the ring around 
the TC from the slightly extincted PMS stars toward the 
TC core, around the most MS massive star 9 1 Ori C. 

In addition, instead of the condition of logN# >22.0 
cm -2 used by Feigelson et al. (|2005|) and Grosso et al. 
(2005), we propose that the condition logN#>22.5 cm" 2 
is more appropriate to isolate the most embedded pop- 
ulation of low-mass PMS sources in the ongoing massive 
star-forming regions in Orion. Fig. 1 shows that the pop- 
ulations of sources with logN# <22.5 cm' 2 (upper and 
middle panels) are distributed throughout the field, includ- 
ing the regions without molecular emission. However, the 
sources with logN# >22.5 cm" 2 (lower panels) are pre- 
dominantly concentrated in the region with CN molecu- 
lar emission, suggesting their physical association with the 
molecular cloud. 



3.2. Stellar densities 

From spatial gridding we can derive the stellar densities 
in terms of the number of stars per cell. However, as ex- 
plained in section 2.1., this method may underestimate the 
values of the density peaks. For that reason, we used the 
close-neighbors method to calculate the densities in the 
low-mass PMS star clusters. In Fig. 3 we report the his- 
tograms of the number of 0.03 pc x 0.03 pc cells centered 
on COUP stars as a function of the number of stars found 
within each cell. The different colors and line types indicate 
the extinction ranges considered in our study. Consistently 
with our results from Fig. 1 and 2 (from spatial gridding), 
Fig. 3 (from the close-neighbors method) shows that the 
cells with the largest number of neighbors correspond to 
the X-ray sources detected toward the TC, the OHC, and 
the OMC1-S region independently of the spatial gridding. 
The densities obtained using cells centered on COUP stars 
are slightly higher than those calculated using the spatial 
gridding, and reflect the real value for the density peaks 
more reliably. 

The surface and volume stellar densities derived for 
the three extinction ranges toward the main density peaks 
(OHC, OMC1-S, the TC envelope, and the TC core) are 
summarized in Table 1. The surface densities were calcu- 
lated by using 0.03 pc x 0.03 pc cells, and the volume den- 
sities using spheres with radius of 0.015 pc, in both cases 
centered on COUP sources. The errors were estimated by 
considering ±0.2 star cell -1 as background/foreground con- 
tamination (section 2.4.). 

The highest stellar density, with nine stars in the 0.03 
pc cell in the range logN# >22.5 cm" 2 , is found in two cells 
located within the OHC; one centered on COUP 621, likely 
a low-mass companion of the infrared source n (Lonsdale 
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Fig. 3. Number of cells versus the number of stars located 
inside the cell, considering 15"xl5" (0.03 pcx0.03 pc) cells 
centered on each COUP X-ray star. Different colors and 
dashed lines show different extinction ranges. The densest 
cells (i.e., cells with the highest number of X-ray stars) are 
located in the OHC, OMC1-S, the TC core and the TC 
envelope (see text). 



et al.. [15521 Grosso et al.. [2lTO5] . and the other on COUP 
622, located only 3.3" south of source n. 

These densities must be considered as lower limits 
since they do not consider the unseen population behind 
very high extinctions. This is particularly critical for the 
extinction range logN# >22.5 cm -2 (see section 2.5.). 
Considering the incompleteness due to obscuration, the sur- 
face and volume stellar densities in the OHC and OMC1-S 
could be as high as (7.0± 0.2)xl0 5 stars pc" 3 and (1.30± 
0.02) xlO 6 for OMC1-S and the OHC, respectively. These 
densities may be even higher due to the presence of unre- 
solved binaries. These values, along with the one found by 
Martin-Pintado et al. ([2005]) in Cep A HW2 (1.6x10 s stars 
pc -3 ), are the highest found in all observed star clusters 
(Bally & Zinnecker, [1003). 

In Fig. 4 we show the derived stellar surface density of 
PMS stars for the three extinction ranges versus the dis- 
tance to the X-ray cluster centers. Left panel shows that 
the TC is dominated by slightly embedded sources concen- 
trated in the central ~ 0.03 pc, where the four MS mas- 
sive stars of the Trapezium are located. The TC envelope 
formed by moderately extincted sources is located at a ra- 
dius of ~ 0.05 pc from the cluster center. The middle and 
right panels show that the embedded populations in the 
OHC and the OMC1-S PMS clusters are very compact with 
sizes ~ 0.03 pc. 

3.2.1. Effect of the cell size 

Table 2 compares the number of sources found in the main 
stellar density peaks (TC envelope, TC core, OMC1-S, and 
OHC) using the close neighbors method with different cells 



We calculated the stellar densities counting the numbers of 
stars contained in rings with widths of ~0.016 pc, centered on 
the cluster centers. 
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Table 1. Number of COUP stars in 15" x 15"(0.03 pcx0.03 pc) cells centered on COUP sources (close-neighbors method), 
and derived stellar densities (surface and volume) for the density peaks seen at different extinctions. 



Region 


Column density range 


# Sources 


Surface stellar density S* 


Volume stellar density p* 




(cm -2 ) 




(stars pc~ 2 ) 


(stars pc — ) 


TC envelope 


<22.0 


5 


(5.5± O.SJIO 3 


(3.5± 0.2)-10 b 


TC core 


22.0-22.5 


i 


(4.4± 0.3)-10 3 


(2.8± 0.2)-10 5 


OMC1-S 


>22.5 


6 


(6.6± 0.3)10 3 


(4.2± 0.2)-10 5 


OMC1-S (corrected) 


>22.5 


10 


(1.10± 0.03)-10 4 


(7.0± 0.2)-10 5 


OHC 


>22.5 


9 


(1.10± 0.03)-10 4 


(7.0± 0.2)-10 5 


OHC (corrected)" 


>22.5 


19 


(2.10± 0.03)10 4 


(1.30± 0.02)-10 6 



a Stellar densities corrected from incompleteness (see section 2.5). 



sizes (0.03 pcx0.03 pc and 0.05 pcx0.05 pc). We defined a 
clustering parameter (a) as the ratio between the number 
of sources found in the 0.03 pcx 0.03 pc cell and the number 
found in the 0.05 pc x 0.05 pc cell. A value of a closer to 
would mean that there is no central concentration of stars 
and hence that the distribution is not clustered. A value of 
a=0.39 would correspond to a uniform stellar distribution. 
A value a=l would mean that all stars are concentrated in 
the smaller cell. Our results show that the clustering pa- 
rameter a is between 0.7 and 0.8 in the clusters of PMS 
low-mass stars found in Orion (Table 2), which means that 
70-80% of the stars in those clusters are concentrated in 
the inner cell of 0.03 pc. If we assume a Gaussian distribu- 
tion for the stellar density p* oc e~( 2r / 7 ) (where r is the 
distance to the center of the cluster and 7 is defined as the 
cluster size), a=0.7-0.8 corresponds to 7 ~ 0.024-0.028 pc. 
This size range is very similar to the size of our smaller 
cell. Consequently, our choice of the cell size (0.03 pc) is 
appropriate for sampling these clusters and for studying its 
stellar densities. 



We note that the cell size is a key parameter for calcu- 
lating stellar densities. The last column in Table 2 shows 
that the volume stellar densities calculated using 0.03 pc 
cells are ~3 times higher than those calculated using larger 
cells of 0.05 pc. This is because the stellar density peaks are 
very concentrated (<0.03 pc). Therefore, the choice of an 
adequate cell size is the key in studying the stellar densities, 
in particular in the very young clusters of PMS low-mass 
stars, which are expected to be very compact. 



3.3. Comparison with other wavelengths 



We compared surveys of the stellar population in Orion 
carried out at different wavelengths with our analysis in 
X-rays. The different catalogs used for the comparison are 
presented in Table 3: one at optical (Hillenbrand, 1997J, and 
three at IR wavelengths (Hillenbrand & Carpenter, 2000 
Muench et al.. 120021 and Lada et al.. l2U04|) . For the compar- 
ison we applied the same spatial gridding as for the X-rays, 
integrating the stellar population in cells of 0.03 pc x 0.03 
pc. Fig. 5 shows the spatial distribution of the stars for the 
different surveys. Note that the color scale representing the 
number of sources per cell is different from that in Figs. 1 
and 2. 



3.3.1. Qualitative comparison with optical and IR 
observations 

We have shown in section 3.1. the importance of considering 
the population of stars in different extinction ranges, to dis- 
tinguish members from the optically visible ONC and from 
the stellar population embedded in the molecular cloud lo- 
cated behind. In this subsection, however, we consider all 
stars observed at different wavelengths (Fig. 5), without 
applying any criteria to separate population at different ex- 
tinctions. This prevents a reliable quantitative stellar den- 
sity analysis study, because we are taking into account stars 
along the same line of sight that belong to different physical 
regions. However, this comparison between different wave- 
lengths can give us interesting additional information. 

The optical observation (upper panel in Fig. 5) reveals 
the ONC population, and presents a morphology consistent 
with the non- extincted X-ray sources (logNjj <22.0 cm~ 2 ; 
upper panel in Fig. 1), with the main density peaks in the 
TC envelope (NE, NW and SE of the Trapezium stars), 
and without a significant population in the TC core. This 
gap in the radial distribution of low-mass stars surrounding 
the most massive star, 9 1 Ori C, was already remarked by 
Hillenbrand (flW]) . 

The three infrared observations show basically the same 
result. Unlike the optical population, the IR population is 
clustered in the TC core, around the most massive of the 
four Trapezium stars (8 1 Ori C) and in the OHC region. 
However, IR surveys do not detect the cluster found in X- 
rays toward OMC1-S. This is probably caused by a very 
high extinction toward this region. 

The clustering in the TC core resembles the morphol- 
ogy presented by slightly extincted COUP sources (22.0 
cm~ 2 <logN#<22.5 cm -2 ; middle panels in Fig 1). This, 
along with the fact that this clustering is observed at IR 
but not at optical wavelengths, suggests that the TC core 
contains more dust than its surroundings in the TC enve- 
lope, consistent with the X-ray results. This might yield 
some clues about the formation mechanism of the massive 
stars in the Trapezium (see the discussion in section 4.2.) 

IR observations toward the OHC region detect both the 
ONC members located along the line of sight and the par- 
tially embedded sources in the molecular cloud, which are 
still visible at these wavelengths. Rivilla et al. (12013b) ) found 
that stars in OHC with logN ff <22.9 cm" 2 (Ay~40 mag) 
are still visible in the IR, while those more embedded are, 
as expected, only visible in X-rays. 

It would be expected that the L' band observations 
(Lada et al.. 12004 1 lower panel in Fig. 5), more sensitive to 
the presence of dust than IR observations at shorter wave- 
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Table 2. Comparison between the density analysis using the close-neighbors method, with cells centered on COUP 
sources with areas of 15" x 15" (0.03 pcx0.03 pc) and 24" x 24" (0.05pcx0.05pc). 



Region 


Column density range 


# Sources 


Clustering parameter 


Volume stellar density ratio 




(cm -2 ) 


15" x 15" 


24" X 24" 


a 


P*[0.03pc] / P*[0.05pc] 






0.03 pc xO.03 pc 


0.05 pcx0.05 pc 






TC envelope 


<22.0 


5 


7 


0.7 


2.9 


TC core 


22.0-22.5 


4 


5 


0.8 


3.3 


OMC1-S 


>22.5 


6 


8 


0.8 


3.1 


OHC 


>22.5 
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13 


0.7 


2.8 
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Fig. 4. Surface stellar density versus the distance to the PMS cluster center for the TC (left), the OHC (middle), and 
OMC1-S (right). The different colors and dashed lines show different extinction ranges: green corresponds to logN# < 
22.0 cm -2 , blue indicates an extinction range of 22.0 cm -2 < logN# < 22.5 cm -2 , and red corresponds to logNjf > 22.5 
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lengths, show a clear clustering in the OHC. However, the 
peak in the OHC region is not as prominent as the peaks 
found in the other two surveys. This is because this region is 
located at the border of the Lada et al. (2004} observation 
(Fig. 5). 

3.3.2. Quantitative comparison with IR observations 

Similar to the procedure used to study the X-ray popula- 
tion in different extinction ranges, Lada et al. (l2004fe used 
the condition J-H<1.5 mag to select partially extincted 
and non-extincted sources (Ay- <8 mag). These authors 
found that the main peak in the TC region (they called 
it Trapezium Core, but it is equivalent to the region con- 
taining our TC envelope and TC core) has a surface stellar 
density of 7.2 x 10 3 stars pc -2 , a value of the same order (al- 
though slightly higher) as the average density of the TC en- 
velope and TC core obtained from the X-ray data (— 5x 10 3 
stars pc -2 ). 

Lada et al. (2004j) also used the condition K-L>1.5 mag 
(roughly corresponding to Ay >15-26, and approximately 
equivalent to our condition logN# >22.5 cm -2 ) to select 
sources deeply embedded in molecular gas and dust. They 
found that the main peak in the embedded population is 
located in the OHC region, with a stellar density of 3.2 x 10 3 
stars pc -2 . This value is a factor 7 below the density we 
have calculated here (2.0xl0 4 stars pc -2 ). This is because 
the most deeply embedded population remains unseen at 
IR wavelengths. Our work confirms the potential of X-ray 
observations to detect the population of heavily extincted 
low-mass stars in massive star cradles. 



3.4. The Trapezium Cluster: intermediate extincted core and 
non-extincted envelope 

Our analysis has revealed for the first time a dense clus- 
ter of X-ray low-mass stars with intermediate extinction 
(22.0 cm -2 <logN# <22.5 cm -2 ) around the most massive 
star in the Trapezium, 6 1 Ori C, surrounded by a non- 
extincted envelope of stars. Most of these stars in the TC 
core were not detected in the optical survey from hillen- 
brand (1997), consistently with their intermediate extinc- 
tion. We have shown that the gap found in the radial dis- 
tribution in the optical survey disappears when considering 
intermediate extincted X-ray stars or IR stars. 

To test whether the COUP stars clustered around 9 1 
Ori C are ONC members with extinctions higher than that 
usually found in the optical visible nebula, or whether they 
are simply slightly extincted members of the background 
OMC cloud, we crossed-checked the COUP catalog with 
optical HST observations of ONC proplyds from Kastner 
et al (|2005p . which can detect stars more obscured than 
the optical Hillenbrand (1997) observations. 

We found that in a close vicinity (20" x 20", -0.04 pc 
x 0.04 pc) of 8 1 Ori C there is an absence of low extincted 
stars, and a cluster of seven X-ray stars with intermedi- 
ate extinction, six of which are proplyds in Kastner et al. 
( 2005) . This clearly confirms that these stars with interme- 
diate extinction are ONC members. 

In Fig. 6 we have plotted the spatial distribution of 
the proplyds with COUP counterparts, distinguishing those 
with logN# <22.0 cm -2 (blue dots in left panel) and 
logN^f >22.0 cm -2 (red dots in right panel). The gray scale 
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Fig. 5. Spatial distribution of stars from four different wavelengths. From upper to lower panels: optical catalog (V and 
Ic bands) from Hillenbrand (|1997|) . IR catalog (H and K bands) from Hillenbrand & Carpenter (|2000p . IR catalog (J, 
H and Kg bands) from Muench et al. (|2002p . and IR catalog (L' band) from Lada et al. p004|) . As in Fig 1., dashed 
contours represent dense gas traced by the integrated intensity emission of CN (N=l-0) (Rodriguez-Franco et al.,[l998). 
The left panels show the position of the X-ray stars as blue dots. Four crosses, an open square and an open triangle show 
the location of the four main-sequence massive Trapezium stars, the OHC, and OMC1-S. The right panels report the 
stellar surface density for the same extinction ranges, derived by counting the number of COUP sources using a kernel 
j^ox of 15" (0.03 pc) (right color scale), superimposed on the CN emission. 
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Fig. 6. Proplyds from the Kastner et al. (|2005p catalog with X-ray counterparts. The gray scale is the 2MASS image at 
K band, and the dashed contours denote the CN N=l-0 emission from the OMC (Rodriguez-Franco et ah. I1998P Left: 
stars with measured values of the hydrogen column density of logN# < 22.0 cm" 2 . Right: stars with measured values of 
the hydrogen column density of logN# > 22.0 cm" 2 . 



is the 2MASS image at K band, where the four massive 
Trapezium stars are clearly visible, and the dashed con- 
tours denote the CN N=l-0 emission from the OMC. Fig. 
6 supports our finding of a core-envelope structure with in- 
termediate extincted stars in the core and non-extincted 
stars in the envelope. 

4. Discussion: Scenarios of massive star formation 
based on accretion theories and comparison 
with Orion observations 

Massive stars (M>8M Q ) evolve very fast and start the nu- 
clear burning when they are still accreting mass from the 
natal molecular core. The radiation pressure from the newly 
formed star has long been proposed to halt accretion, pre- 
venting the formation of the most massive stars (Wolfire & 
Cassinelli,[TWI). 

One of the first theories proposed to resolve this prob- 
lem was the coalescence of low- and intermediate-mass stars 
in a dense stellar cluster (Bonnell et ah. HMH l2"0TJ2l l2"0tI5|) , 
a violent process that does not require subsequent mass 
accretion. The observational evidence suggesting that a vi- 
olent event occurred in the OHC region around 500 yr ago 
(Zapata et al.. 120091 Rodriguez et al, Gomez et al., 2008 
Goddi et ah. 12011)1 . along with the high stellar density re- 
ported in our work (~ 10 6 stars pc~ 3 ) would make coa- 
lescence possible in this region. In Appendix A we discuss 
the possibility of a stellar collision in the high dense clus- 
ter revealed by Chandra in the OHC region, and conclude 
that a stellar interaction could have occurred in the cluster 
although it is uncertain whether this event has finally led 
to the formation of a massive star. 

Stellar coalescence has always been seen as an exotic 
mechanism. This is because it requires higher stellar densi- 
ties (10 6 -10 9 stars pc -3 ) than those usually found in stellar 
clusters to produce enough stellar collisions in 10 5 -10 6 yr to 



form massive stars. Furthermore, this theory lost momen- 
tum when several works based on accretion showed that 
it was possible to overcome the radiation pressure prob- 
lem without invoking coalescence (Yorke & Sonhalter, 2002; 
Cesaroni et al.. l2006p . 

A general picture for massive star formation should 
likely be based on the two theories based on accretion: 
monolithic core accretion and competitive accretion. The 
main difference between these two theories is how the stars 
can accrete the mass needed to become massive. 

According to the monolithic core accretion theory, ini- 
tially proposed by McKee & Tan ((20021 [2003)) . massive stars 
form from massive cores supported against self-gravity by 
its internal turbulence, with accretion rates high enough 
to overcome the radiation pressure. The massive core col- 
lapses monolithically and produces a single massive object 
or a few massive stars, rather than many low-mass stars 
(Krumholz et al.. [20061 [2009)1 . 

The competitive accretion theory relies on the fact that 
nearly all massive stars are formed in stellar clusters with 
low-mass stars (Clarke, Bonnell & Hillenbrand, 120001 Lada 
& Lada, 2003)). This theory assumes that the parental 
molecular cloud fragments into several condensations with 
masses around the Jeans masa3. Numerical simulations 
have shown that gravitational fragmentation of a molecular 
cloud produces condensations with masses of about Mj eans 
< 1 Mq (Klessen et al.. [T998l Bate et ah. [20041 Bonnell et 
al., I2004p that collapse and form a low-mass star cluster. 
An additional mechanism is then needed to explain the for- 
mation of massive stars. It has been proposed that this 
mechanism is continuous gas accretion. Because fragmen- 
tation is inefficient in accumulating gas onto stars, there is 
a common gas reservoir in the core that can potentially be 
accreted (Bonnell et al.. 12003)) . Observations of pre-stellar 

2 The Jeans mass, Mj ean s, is the mass necessary for an object 
to be gravitationally bound against its thermal support. 
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Table 3. Surveys of Orion at different wavelengths. 



Wavelength 


Telescope 


Reference 


optical (bands V and Ic) 


KPNO 0.9 m + data from previous literature 


Hillenbrand (1997) 


IR (H and K bands) 


Keck I 10m 


Hillenbrand & Carpenter (l2000h 


IR (J, H and K 3 bands) 


1.2 m telescope at FLWO & 3.5 m NTT in La Silla 


Muench et al. (I2002D 


IR (L') 


VLT UT1 


Lada et al. (1200411 



structures and young clusters indeed show that most of the 
total mass (Krumholz & Bonnell, 2007) is in a distributed 
gaseous form, so that it can potentially be accreted. Then, 
after low-mass stars are formed, they fall into local gravita- 
tional potential wells forming small-N clusters. This overall 
cluster potential well funnels gas down to the potential cen- 
ter, where it is captured via Bondi-Hoyle accretion by the 
"privileged" stars located at the potential center, progres- 
sively increasing their masses. 

Clearly, core accretion and competitive accretion are 
rather different ways to explain massive star formation. 
However, it is not unreasonable to think of a scenario that 
involves a combination of the two theories, whenever the 
physical conditions allow them to work. Some authors have 
already suggested mixed scenarios for massive star forma- 
tion (Peretto et al.. [2"DTJ^1 Wang et al. l2TJTU| . 

In this section we present three possible scenarios for 
massive star formation, where in principle both accretion 
theories could work. In practice, we show that the domi- 
nant mode of star formation is determined by the level of 
fragmentation in the core. While some works claim that 
radiation feedback from the stars can stop fragmentation 
(Krumholz et al., 2006), and allows monolithic collapse to 
work, others claim that the core inevitably fragments into 
smaller condensations to form many stars (Clark & Bonnell, 
[211061 Dobbs et al., [Il l Bonn ell & Bate,[Hl Peters et al., 
120101 Federrath et al., I2010[) . In this context, our scenar- 
ios consider different levels of fragmentation: i) low frag- 
mentation and monolithic core accretion; ii) intermediate 
fragmentation and subcore accretion; iii) high fragmenta- 
tion, subcore accretion and subsequent competitive accre- 
tion. We compare the predictions from each scenario with 
the results of our analysis presented in section 3. 

4.1. Low fragmentation: Monolithic core accretion 

In this scenario the radiative feedback halts fragmentation 
of the core, which collapses monolithically, allowing the for- 
mation of a single massive (or a few) objects. This is the 
core-to-star picture proposed by Krumholz (|2006l I2009|> . 
This general collapse would rapidly consume the available 
gas, so that although subsequent Bondi-Hoyle accretion 
of the remanent gas could occur once the massive star is 
formed, the final mass is mainly determined by the core- 
accretion process. A simple schematic picture of this sce- 
nario is presented in the upper row in Fig. 7. 

According to this scenario, massive stars should be born 
in isolation or be accompanied by other massive stars, 
rather than with a low-mass stellar cluster inside the core. 
Therefore, although the theory could explain the small frac- 
tion of massive stars found in isolation (de Wit et al., 2004), 
it would not explain the presence of low-mass stellar clus- 
ters with the massive star cradles in Orion (TC, the OHC 
and and OMC1-S; see section 3), and also in most stel- 



lar clusters (Clarke, Bonnell & Hillenbrand, 2000; Lada & 
Lada. [20031) . 

In addition, it is known that massive cores in infrared 
dark clouds (IRDCs), believed to represent the initial con- 
ditions for massive stars and star cluster formation, are 
usually found to fragment into smaller condensations of 
low- to intermediate masses (see e.g. Wang et al., I201ip . 
which is in contrast with the idea of a massive star being 
formed from a single massive core. Other works (Bonnell & 
Bate. |2TjT)51 Dobbs et al.. I2005[) have suggested that such a 
massive core, even if supported by turbulence, very likely 
fragments to form many stars. 

4.2. Intermediate fragmentation: subcore accretion 

One could think of an scenario based on monolithic col- 
lapse, but acting at subcore scales. In this scenario, the 
core would suffer fragmentation, generating several sub- 
cores with different masses and sizes. The subcores with 
lower masses would produce low-mass stars, while those 
with higher masses could potentially form massive stars, in 
both cases via subcore accretion. A schematic picture ex- 
plaining this scenario is presented in the middle row in Fig. 
7. In principle, massive stars could be formed in isolation 
or be accompanied by low-mass stars, depending on how 
the core initially fragments into subcores. 

Once the stellar cluster is born, subsequent Bondi-Hoyle 
accretion from the gas reservoir still available in the core 
would be possible. However, it is difficult for competitive 
accretion to work in this context because there is no overall 
well-defined central cluster potential well, but a collection 
of local potential wells around each formed star. Therefore, 
there is no efficient funneling of gas to the center of the 
stellar cluster. 

This scenario is very similar to the results presented by 
the simulations by Krumholz et al. (|20111 12012[) and it is 
able to produce stars with different masses. These works 
implement the theory of core accretion in the general col- 
lapse of a cloud instead of the single massive core considered 
in previous simulations (Krumholz et al., 2006, 2009). The 
outcome of their simulations includes some massive stars 
formed nearly in isolation (consistently with their previous 
works) and also others located in low-mass stellar clusters, 
which would be more compatible with our results in Orion. 

However, this scenario presents several problems. 
Krumholz et al. (|2011l [2012) argued that radiative feedback 
from the stars stops further fragmentation, explaining why 
not all subcores reach low-masses ~ 1 M Q during the col- 
lapse phase. However, it is not clear that radiative feedback 
from the newly stars could be so efficient to stop fragmen- 
tation. Several other works point in the opposite direction, 
i.e., massive cores should fragment into many stars (Dobbs 
et al.. 120051 Clar k fc B onnell, [20061 Bonne ll & Bate. [20061 
Federrath et al.. 120101 Peters et al.. l2010p . especially with 
the presence of outflows that are not included in the simu- 
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Fig. 7. Scenarios for massive star formation based on the two accretion theories. The dominant role of each theory is 
mainly determined by the level of fragmentation in the core (see text). Upper panels: low fragmentation and mono- 
lithic core accretion. Middle panel: intermediate fragmentation and subsequent subcore accretion. Lower panel: high 
fragmentation, subcore accretion forming stellar seeds, and competitive accretion. 



lation by Krumholz et al. (20111 l2"t)12p and could favor the 
fragmentation (Cunningham et al.. I2011[) . 

In addition, there is other evidence indicating that this 
scenario does not fully explain, at least, the massive star for- 
mation that occurred in Orion. In this scenario fragmenta- 
tion does not always produce more massive subcores at the 
center of the initial core (see Krumholz et al. 120111 [5012). 
Therefore, one would not expect a trend for massive stars 
to form at the center of the core, as seen in the TC (see 
discussion in section 4.3). 

4.3. High fragmentation: subcore accretion forming stellar 
seeds and subsequent competitive accretion 

We assume now that the core fragments into subcores, 
which collapse producing a cluster of low-mass stars via 
subcore accretion. Only a fraction of the gas is directly 
incorporated into the low-mass stars, while the rest is dis- 
tributed throughout the core. In this context, the stars com- 
pete to accrete this material. The stars located near the cen- 
ter of the potential well created by the whole cluster benefit 
from the higher gravitational attraction and gather matter 
via Bondi-Hoyle accretion at much higher rates, becoming 
higher mass stars. The stars that are not at the center of 
the potential well do not accrete significant amounts of gas, 
and remain low- and intermediate-mass stars. This scenario 
with high fragmentation (lower panel in Fig. 7) is essen- 
tially coincident with competitive accretion, and naturally 
produces a whole cluster of stars with different masses. 



Bonnell et al. (|2004p simulated the evolution of a tur- 
bulent molecular cloud with properties similar to those in 
Orion, where fragmentation and competitive accretion play 
the main role in the formation of massive stars. The result 
of their simulation (see their Figure 5 and also Fig. 14 in 
Zinnecker & Yorke. [2007[) shows a ~ 0.6 pc x 0.6 pc region 
where the massive stars are formed in the center of low- 
mass star clusters. The similarity between the morphology 
of the whole cluster with the Orion observations presented 
in this paper is remarkable (Fig. 1). Moreover, there is a 
trend in the TC for higher mass stars to be found at closer 
distances to the most massive star 9 1 Ori C. Fig. 8 shows 
the mass of the stars in the TC versus the distance to 9 1 
Ori C (data from the optical observations by Hillenbrand, 
1997). This mass segregation is primordial, i.e., the mas- 
sive stars were initially born in the cluster center without 
migration (Bonnell, Larson & Zinnecker, |2QQ7l Reggiani et 
al. . 1201 1|) . Therefore, this suggests a preference for the most 
massive stars to be born in the densest and central parts 
of the low-mass star cluster, as expected in the competitive 
accretion theory. 

This scenario can also naturally explains the structure of 
extincted and non-extincted low-mass stars found in the TC 
(Fig. 1, 2 and 6). The non-extincted stars in the envelope of 
the TC could have lost the battle of competitive accretion - 
keeping their low-mass - against the four Trapezium stars in 
the core, which would have grown until they reached their 
current high masses. The non-extincted PMS stars in the 
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envelope of the TC are therefore devoid of gas while those 
in the core are still partially embedded in gas and dust. 

Given that this scenario provides the best explanation 
for our results, we evaluated whether Bondi-Hoyle accretion 
can indeed gather enough mass to produce a massive star 
in the Orion environment. The accretion rate m is (Bonnell 
& Bate, [5DM1) 



, (GM) 2 
m = 4np — . 

1 1*5 



(1) 



where p is the volume gas density, G the gravitational 
constant, and v the relative velocity between the gas and 
the accreting star. Considering that the relative gas velocity 
follows the turbulent velocity size-scale relation v ~ R 1 / 2 
(Larson et al., 119811 Heyer & Brunt, 2004), and assum- 
ing that the gas density p is centrally condensed with 
p r~j i?" 1 (as found in other massive star-forming regions 
like Monoceros R2, Choi et al., [2000), the relation be- 
tween the velocity and the mass of the core is v ~ M 1 / 4 . 
Combining this with eq. 1 and integrating this expression, 
one can obtain the time needed to build up a massive star 



t = t. 



1 



1 



irpG 2 M? /4 Im? 



/4 Afl/4 



(2) 



where M is the mass of the star and Mj is the initial 
mass of the "stellar seed", and t see d is the time needed to 
form the stellar seed with initial mass M,. 

The relative velocity of the gas v is a key parameter 
to determine whether competitive accretion is efficient or 
not (Krumholz et al. . [211051 Bonnell & Bate. I2006"|) . It has 
been suggested that turbulence injected by winds, proto- 
stellar outflows, or ionized regions could increase the rel- 
ative velocity of the gas, preventing competitive accretion 
from being the dominant process in the massive star forma- 
tion (Krumholz et al.,[2005). However, the results presented 
in this work have shown that massive stars in Orion have 
formed in small-N clusters (tens of stars, see Table 1), where 
the gas velocity dispersion is expected to be lower (Bonnell 
et al. 12003112005) . 

Assuming a Maxwellian distribution of velocities, we 
used the root mean square of the gas velocity, vrmSi as a 
good estimate of the relative velocity v. The value of vrms 
is related to the linewidth of the molecular lines, Av, with 
the expression 



vrms 




2Ln2 



Av, 



(3) 



We use here the linewidth of vibrationally excited HC3N 
(HC3N*, hereafter) because this molecule is particularly 
well suited to study the kinematics of the dense gas close to 
very young massive stars (Wyrowski et al. . 119991 de Vicente 
et a l [20001 [2002]) . The OHC observa tions by de Vicente et 
al. (|2002) and Rivilla et al. (|2013cl in preparation) show 
that Av ~ 7 km s" 1 . This value is consistent with those 
observed with other dense molecular tracers as NH3 (~5 
km - 1 , Go ddi et al.. ISOLD) or CH 3 OH (4-6 km s -1 , Peng 
et al. I2012|) . Using eq. 3, the HC3N* linewidth is equivalent 
to VijA/s=5.2 km s . 

Fig. 9 shows the time needed to form a massive star 
from stellar seeds with M^=l M Q by subcore accretion and 




Fig. 8. Mass of the stars in the TC versus the distance 
to the most massive star 9 1 Ori C (data from the optical 
observation presented by Hillenbrand. 11997]) . 



subsequent competitive accretion. The blue line indicates 
the accretion time from the expressions from Krumholz et 
al., ([20L2)) and McKee & Tan ([20031) to form the stellar seed. 
The red line corresponds to the Bondy-Hoyle accretion time 
using eq. 2, with v = urms=5.2 km s . The volume- and 
surface densities used are p=1.95x 10" 16 g cm" 3 and £=12 
g cm -2 , calculated from the average gas density in the OHC 
of N[H 2 ]=5xl0 7 cm" 3 (Morris et al. fT980]> . According to 
Fig. 9, a 10 Mq star (approximately the mass estimated for 
the massive objects found in this region) can be formed in 
the OHC by accretion of gas in < 4x 10 5 years. We note that 
this accretion time estimate does not demonstrate by itself 
that massive stars are formed by competitive accretion in 
the OHC, but it shows that it could be possible. 

In the OMC1-S region, where the gas density is lower 
but on the order of the OHC region (see CN contours in Fig. 
1 and 2), and where a similar (or lower) turbulence may be 
expected, massive stars could form in the same way. 

In the case of the TC, the current gas density is ^10" 17 
g cm" 3 (Hillenbrand & Hartmann, 1998 ; Bonnell & Bate, 
2006), one order of magnitude lower than that in the OHC 
(~1.95xl0~ 16 g cm" 3 ). This implies that competitive ac- 
cretion could be less important according to eq. 1. However, 
as indicated by Krumholz et al. (|201ip . this density must 
have been higher by an order of magnitude at the moment 
where the stars were formed, because the cluster has likely 
expanded due to gas ejection (Kroupa et al., 120011 Tan 
et al.. 12006]) . Then, competitive accretion could still have 
been responsible for the formation of the massive stars in 
the Trapezium. 

In summary, from our combined data (CN molecular 
emission and X-ray PMS stars), it seems that at least two 
basic ingredients are needed to form massive stars: i) dense 
gas, and ii) a cluster of low-mass stars. It is remarkable 
that the low-mass star clusters appear only at certain posi- 
tions of the north-south elongated ridge of dense molecular 
material (Fig. 1), suggesting that dense molecular mate- 
rial only does not necessarily form massive stars. This idea 



12 



V. M. Rivilla et al.: The role of low-mass star clusters in massive star formation. The Orion Case. 




10 1 1 10 1 10 2 

M (M ) 

Fig. 9. Accretion time necessary to form a massive star in 
the OHC in the high fragmentation scenario. The blue line 
corresponds to the time needed to form the stellar seeds of 
1 M Q via subcore accretion. We have used the expressions 
from Krumholz et al. p012p for the accretion rate. The 
red line corresponds to the accretion time needed to reach 
the final mass M, using eq. (2). The volume- and surface 
densities used are p—1.95x 10~ 16 g cm -3 and £=12 g cm~ 2 . 



has also been suggested by Krumholz et al. (2012). Along 
with dense gas, our results suggest that the presence of 
low-mass star clusters is also a necessary but not sufficient 
condition for massive star formation. This is supported by 
the existence of dense low-mass star clusters, like the one 
at ^[83.84, -5.395] (see Fig. 1), which does not present any 
massive star. This low-mass star cluster is seen in the non- 
extincted column density range (logN H <22.0 cm" 2 ) in the 
X-ray observations (upper panel in Fig. 1), and also in the 
IR surveys and in the optical observation by Hillenbrand 
(| 19971 see upper panel in Fig. 5), confirming that it belongs 
to the ONC. In this region, the lack of dense molecular gas 
suggests that this low-mass star cluster was unable to ac- 
crete additional material to increase its mass to form more 
massive objects. 



5. Conclusions 

We have analyzed the distribution of PMS low-mass stars 
with X-ray emission in Orion observed by Chandra as a 
function of extinction, and calculated their stellar densities. 
Our results show that PMS low-mass stars cluster toward 
the three regions connected to massive stars: the Trapezium 
Cluster (TC), the Orion Hot Core (OHC), and the OMC1- 
S region. We derived PMS low-mass stellar densities of 10 5 
stars pc" 3 in the TC and OMC1-S, and 10 6 stars pc -3 in 
the OHC. These very high stellar densities and the close 
association between the low-mass star clusters with the re- 
gions where massive star formation has recently occurred 
(TC) or it is currently taking place (OHC and OMC1-S) 
suggest that not only dense gas is needed to form massive 
objects, but also a cluster of low-mass stars. 



The X-ray observations show for the first time in the 
TC that low-mass stars with intermediate extinction are 
clustered toward the position of the most massive star 6 1 
Ori C, which is surrounded by a ring of non-extincted PMS 
low-mass stars. This 'envelope-core' structure is also sup- 
ported from infrared and optical observations. 

We proposed several scenarios for massive star forma- 
tion in which the theories based on accretion (core accre- 
tion and competitive accretion) play different roles and the 
parental core presents different levels of fragmentation: i) 
low fragmentation and monolithic core accretion; ii) inter- 
mediate fragmentation and subcore accretion; iii) high frag- 
mentation, subcore accretion and subsequent competitive 
accretion. We compared the predictions of each scenario 
with the findings of our analysis of the X-ray observations 
by COUP. 

The first scenario, based on the classic core accretion 
theory (core-to-star), does not explain the formation of a 
significant population of low-mass stars, contradicting the 
results presented in this paper. 

We showed that the second scenario (intermediate frag- 
mentation) is not fully appropriate to explain our results. 
In this scenario fragmentation does not always produce 
more massive subcores at the center of the parental core. 
Therefore, we would not expect the trend for massive stars 
to form at the center of the core that we observe in the 
TC. Moreover, the presence of a significant low-mass stellar 
population revealed by X-rays in the massive star cradles 
suggests that the level of fragmentation of the initial core 
is high. 

This led us to a third scenario, with higher fragmen- 
tation forming a low-mass stellar cluster via subcore ac- 
cretion. These "stellar seeds" can gather enough mass by 
competitive accretion, forming more massive stars at the 
cluster center. Therefore, the outcome of this scenario pre- 
dicts our findings. Furthermore, we showed that the forma- 
tion of massive stars in a core with the physical conditions 
measured in the OHC can be explained in terms of com- 
petitive accretion. Competitive accretion could also work in 
OMC1-S, with similar gas properties, and in the TC, which 
was denser at the moment when the stars were formed. 

Finally, given that the stellar density in the OHC is 
one of the highest reported in the galaxy (1.30xl0 6 stars 
pc~ 3 ) and because there is much observational evidence 
that a violent process occurred in the core, we discuss the 
possibility of a coalescence event in the OHC stellar cluster 
in Appendix A. We conclude that although coalescence does 
not seem a common mechanism for building up massive 
stars, a single stellar merger may occur in the evolution of 
the cluster of the OHC, favored for the presence of disks, 
binaries, and accretion. 
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Appendix A: Stellar collisions in dense clusters 

The coalescence of low- and intermediate-mass stars has 
been proposed as one of the mechanisms that may 
form massive stars (Bonnell et al., 119981 Stahler et al., 
l2DUm Bonnell & Bate, [20021 Bally & Zinnecker, l2TjU5)) . 
Observations have shown that stars are almost never born 
in isolation but favor formation in clusters (Lada & Lada 
2003). Bonnell et al. (|T9"9"8"1) and Stahler et al. (j2TjDTj|) ar- 
gued that high stellar densities may exist for a brief period 
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A. 1.2. Star-(star+disk) encounter 

In this section, we extend the previous case by adding the 
presence of a circumstcllar disk. Calculating the gravita- 
tional potential generated by a disk is difficult in general 
(Binney & Tremaine, 1987 ') . For this reason, we made a rea- 
sonable approximation to characterize the disk, so that the 
calculation of its potential is rather simple. We consider the 
disk to be an infinitesimal thin surface, with surface density 

*)-» -. ,; + V (A - 2) 

2tt (r 2 + r 2 d ) 

where /3 = (1 - l/v^)" 1 , r d and M d are the radius and 
mass of the disk, and r is the radial distance to the disk 
center. Applying the Gauss theorem, one can obtain the 
gravitational potential created by the disk as 

v r + r i 

This potential generated by the disk adds one term to 
the expression of the stellar collision rate t~ ^ : 



of time during the very embedded early evolution of stellar 
clusters in molecular cores, making stellar collisions more 
likely. Therefore, interactions could occurr among members 
of a protostellar cluster, especially at the early stages of 
cluster evolution. These interactions could play a role in 
the formation of high-mass stars. 

Coalescence requires stellar densities in the order of 10 6 - 
10 8 s tars pc~ 3 (Bonnell & Bate, [20051 Bally & Zinnecker, 
2005) that are higher than the highest densities (of ~ 10 5 
stars pc~ 3 ) found observationally (e.g. Claussen et al. . 119941 
or Figer et al., I2002[) . However, we have detected in this 
work a very high stellar density in the OHC region of 10 6 
stars pc -3 , a range where stellar collisions may occur. This 
indicates that a collision could be possible in the OHC. This 
idea is also supported by some observational evidence that 
suggests that a violent event occurred in the OHC in the 
past: a large-scale CO outflow (Kwan & Scoville, 1976), the 
H 2 fingers (Taylor et al.. 1155^ Al len & Burton, 11993), the 
CO filaments (Zapata et al.. l2009|) . and the proper motions 
measured for the massive objects BN, I and n (Rodriguez 
et al, I2TI051 Gomez et al., [2T7TM Goddi et al., l2"0TTj) . all 
of them pointing toward approximately the same central 
position within the OHC region (see Fig. 1 of Zapata et 
al.. l2UTl"j) . Rodriguez et al. ([M)5) and Zapata et al. ([21309")) 
suggested that this violent event could have been triggered 
by an interaction among several stars in the protostellar 
cluster. This interaction is consistent with a close encounter 
that would induce stars to coalescence in this particular 
high-density cluster revealed by X-ray observations. 

In this appendix, we derive some expressions to describe 
stellar collisions in a cluster and apply them to the prop- 
erties of the OHC with the aim to determine if a stellar 
collision was possible in the region. We also compare our 
results with other works that include dynamical simulations 
of stellar encounters. 

A.l. Physics of stellar collisions 

We follow Binney & Tremaine 1 1 98 71 to study the case where 
two naked stars encounter (A. 1.1) and extend this case to 
a encounter between a star and a star with a circumstellar 
disk (A. 1.2). In section A. 1.3 we present the possibility of 
disk-induced stellar capture in an encounter forming a bi- 
nary, and finally in A. 1.4, we give a simple analysis of the 
orbital decay of a binary trigger by accretion onto the sys- 
tem, which may produce a collision between the two stars. 

A. 1.1. Direct star-star encounter 

We consider an encounter between two stars with mass M. 
Applying conservation of momentum and energy in the sys- 
tem, it is possible to obtain the collision rate (essentially 
the inverse of the time needed to produce a collision, t co u) 

-i- [yr -i]=4V^r c 2 oH (l + ^), (A.l) 

'coll ^ * / 

where p„ is the volume stellar density and c* is the 
velocity dispersion of the cluster. The first term corresponds 
to the geometrical cross section of the stars, and the second 
is due to gravitational focusing. The two stars collide when 
the collision distance r co u is equal to the distance between 
the centers of mass of the stars, i.e., the sum of their radii: 
r co z; = Ri + R 2 - 



1 r -n , r- 2 (, GA/ f3GM d \ . . 

— [yr !] = A^p*ay coll 1 + ^— + P (A.4) 



In this case, a collision in the system would occur when 
Tcoii — rdisfc, considerably increasing the probability of col- 
lision. 



A. 1.3. Capture of the incident star by the system star+disk 

When a collision as described in A. 2. is produced, the disk 
may capture the incident star. The condition for a star to be 
captured by the disk is that the energy of the incident star 
is lower than the potential energy created by the star+disk 
system: 

1 2 PGM d m GMm /A rN 

Z V r coll+ r d TcoU 

This gives a maximum value of t co u for a star to be 
disk-captured in the collision: 

W = -z + Mj . (A.6) 

Then, if t co ii < r raM , a binary could be formed as a 
consequence of the collision. 



A. 1.4. Orbital decay in the binary system induced by 
accretion 

The subsequent accretion onto a binary system will harden 
the binary, moving the orbit closer. This produces an orbital 
decay (Bate et al., I2002j) and eventually could end in a 
collision between the two stars. The evolution of the semi- 
mayor axis of the binary orbit (a) due to accretion can be 
described as (Rayburn, I1976[) 

d=-3a^, (A.7) 
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where we considered that the two binaries have the same 
mass M. Integrating this expression from the moment when 
both stars have initial masses M, and an initial semi-major 
orbit a,, to the moment when the stars have accreted mass 
to form stars with mass M/, we obtain the value of the final 
semi-major axis 

af = ai {w) <ai - (A - 8) 

A.2. The OHC stellar cluster 

In section A.l. we have derived the expressions for the time 
needed for a stellar collision to occur in a stellar cluster, 
t C oii- Fig. A.l. shows the collision time t co u as a function 
of the stellar velocity dispersion <7* in a cluster of low-mass 
stars with a stellar density equal to that found in the OHC 
of 1.30xl0 6 stars pc~ 3 . We considered collisions between 
1 M stars (left panel), and 10 M stars (right panel). 
We included the case of a direct collision between 'naked' 
stars (eq. A.l.) and the case with the likely presence of 
circumstellar disks (eq. A. 2.) (Eisner et al., 2005}. In the 
direct 'naked' star collision, the collision distance r co « is 
the sum of the radii of the stars. We considered that the 
radii of the stars are those of stars with 1 M and 10 M 
at 10 5 yr (left and right panels, respectively), the typical 
timescale for massive star formation. In the case including 
disks, the collision distance can be as high as the size of the 
disk (Bally & Zinnecker, 2005). We considered disks with 
a mass of 0.1 times the mass of the star and reasonable 
radii of 10 and 100 AU, which can significantly reduce the 
collision time via disk-captures (Zinnecker & Yorkc, 2007 
Davies et al. . 120061 Bonnell & Bate. 12005]) . Moreover, bina- 
ries and/or multiple stellar systems are commonly found in 
Orion (Kraus et al.. 12009} and could additionally decrease 
t co u, favoring coalescence. 

There are no direct measurements of the stellar veloc- 
ity dispersion of the OHC cluster ov However, lower and 
upper limits can be estimated. As indicated by Bonnell & 
Bate (I2006p . one might expect that the stars forming part 
of a small-N young cluster like that in the OHC have a 
low-velocity dispersion. These authors proposed to use val- 
ues as low as 0.4 km s _1 . Gomez et al. (2008) claimed that 
typical random motions of recently formed stars have a ve- 
locity dispersion of 1-2 km s _1 . The ONC has a measured 
velocity dispersion of ^2.3 km s _1 (van Altena et al.,[l988). 
On the other hand, we can use the gas velocity dispersion 
of the OHC as a conservative upper limit to the stellar ve- 
locity dispersion, because stars are expected to be much 
less affected by turbulence than gas. The gas velocity dis- 
persion in the OHC, using HC3N* (see section 4) is a = 
Av/(2V2 Im'2) ~3kms^ (de Vicente et al. . 120021 Rivilla 
et al.. l2013cl in preparation). This upper limit for cr* is very 
consistent with the values presented before. 

According to Fig. A.l., a stellar cluster as dense as the 
OHC with a stellar velocity dispersion of c* ~ 3 km s _1 
could experience one stellar collision of two 1 M stars with 
circumstellar disks with radius 10-100 AU (left panel) in a 
collision time of about ~10 5 yr, making it plausible that the 
OHC cluster has suffered an encounter of low-mass stars. If 
a 10 M star (which can be formed by competitive accretion 
in the low-mass star cluster, see section 4.3) is involved, 
the collision time is even shorter. Therefore, this analysis 



suggests that the violent event that has been proposed for 
the OHC could be a collision between two members of the 
dense cluster revealed by X-rays studied in our work. 

At this point, one may wonder if this close encounter 
may lead to a direct star-star merger. This is a necessary 
condition (although it may not be sufficient) to build up 
a more massive star. Of course, a rigorous characterization 
of the evolution of this system would need complex simula- 
tions that are out of the scope of this appendix. However, 
we can make some assumptions and simple calculations to 
determine if this is, at least, viable. 

After one star has impacted the disk of a neighboring 
star, it can be captured to form a binary if its energy is 
lower than the potential energy created by the star+disk 
system (Bally & Zinnecker. 120051 see also section A. 1.3.). 
This condition sets a maximum collision distance at which 
the star can be captured, r max . According to eq. A. 6. this 
maximum distance in a cluster with stellar velocity disper- 
sion of 3 km s _1 and disk mass Md=0.1M star , is ~250 AU 
for a 10 M Q star and -2500 AU for a 1 M star. Therefore, 
the encounters we considered in this work with r co u <100 
AU can potentially form a binary system as a consequence 
of the disk capture. 

Once the binary system is formed, subsequent accretion 
of the parental molecular core will harden the binary, trig- 
gering an orbital decay (Bate et al.. 12002} . Bonnell & Bate 
(2005) proposed that competitive accretion in a dense clus- 
ter naturally forms more massive and closer binaries. If we 
consider stars with 1 M with disks of reasonable radii of 
—100 AU, collisions with r co jz=r disk are expected to occur 
in 7x 10 4 yr for a stellar velocity dispersion of ~ 3 km s _1 
(Fig. A.l.). Fig. 9 shows that the binary system, if located 
at the center of the OHC low-mass star cluster, can gather 
a mass of 5 M via competitive accretion in ~2-3xl0 5 yr. 
Therefore, in several 10 5 yr the final individual masses of 
the binary would be 6 M . Then, from eq. A. 8. we find that 
the continuous accretion onto the binary could reduce the 
semi-major axis of the orbit in a factor 0.005. Considering 
an initial semi-major axis ai = 100-1000 AU, &f would then 
be 0.5-5 AU. These simple calculations are roughly consis- 
tent with the results from numerical simulations of Bonnell 
& Bate (2005; sec their Fig. 2), where the binaries forming 
in this way typically have very eccentric orbits with peri- 
astron separations comparable to the stellar radii, allowing 
mergers in clusters with stellar densities of 10 6 stars pc~ 3 . 
The stellar density in the OHC is similar to this value, and 
the final semi-axis binary separation we estimated can be 
as low as tens of solar radii, which is slightly higher but 
on the order of the sum of the radius of two 6 M stars 
with ~10 5 yr (Siess et al- 12000} . This implies that a col- 
lision due to orbital decay in a binary system induced by 
accretion could occur in the OHC. 

Whether such a merger leads to a more massive object 
is difficult to asses. However, we can compare our analysis 
with other works where stellar coalescence play a role in 
massive star formation (Davies et al.. 120061 and Moeckel & 
Clarke. [201 ip . The first work suggests that very high stellar 
densities (10 9 M pc 3 ) are needed to form a 50 M star 
through mergers in the adequate timescale. This is three 
orders of magnitude higher than what we found in OHC. 
However, as Davies et al. (120061) recognized in their work, 
this stellar density can decrease significantly if circumstellar 
disks or primordial binaries are considered. Indeed, Davies 
et al. (|2006p showed that including circumstellar disks de- 
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Fig. A.l. Stellar collision time t co u versus stellar velocity dispersion cr*. We have used the lower limit for the stellar 
density we have determined for the OHC of 1.30xl0 6 stars pc -3 . We have considered collisions between 1 M Q stars (left 
panel), and 10 M Q stars (right panel). We have included the case of a direct collision between 'naked' stars (eq. A.l.) 
and the case with the likely presence of circumstellar disks (eq. A. 2) with radii 10 and 100 AU (see text). 



creases the collision time (and the stellar density needed) 
by a factor up of 3-10, and that the presence of primor- 
dial binaries additionally decreases the stellar density by a 
factor of 3. In addition, these authors considered the forma- 
tion of a 50 M Q -star by coalescence, while the most massive 
objects in the OHC region have masses of only ~10 M Q . 
According to Fig. 5 in Davies et al. (2006|). a star with 10 
M Q is formed in a time ~2.5 shorter than a star with 50 M Q . 
Then, the stellar density presented in Davies et al. (2006) 
could be reduced by nearly 2 orders of magnitude, resulting 
in 10 7 M Q pc 3 . This density is still higher than what we 
found in the OHC, suggesting that coalescence cannot be 
considered a general process for building up massive stars. 
This is because to be effective it needs multiple stellar col- 
lisions (Davies et al., I2006p . According to our analysis, a 
single merger could occur in the cluster evolution, but not 
many. 

This also agrees very well with the recent work by 
Moeckel & Clarke pOTT|) . They simulated an O rion-like 
cluster, assuming a stellar density of 2.5 xlO 4 stars pc~ 3 , 
and followed its evolution until 3 Myr, the approximate age 
of the cluster. Even with this relatively low stellar density, 
in two of their five runs they observed two stellar colli- 
sions involving binaries at the central and dense regions 
of the cluster. However, they concluded that an Orion-like 
cluster can be explained without many collisions. This is 
evidence against the need for collision to explain massive 
star formation. Furthermore, even when a collision is in- 
volved, Moeckel & Clarke (2011) found that the final mass 
is mainly determined by accretion growth. Bonnell & Bate 
(2005) also showed that in collisions the mass accreted by 
accretion is generally higher than the mass directly gath- 
ered by the mergers. 



In conclusion, coalescence certainly does not seem a 
common mechanism for building up massive stars. However, 
this does not directly rule out the possibility that it may 
work under certain circumstances, or that a stellar collision 
may occur in a massive star cradle, especially in regions 
with high stellar densities. In particular, there is enough 
evidence to make a stellar collision possible in the OHC. 
This is supported by the very high stellar density we found 
in our work (~10 6 stars pc~ 3 ). Whether this collision pro- 
duced a massive star as a result of a stellar merger still 
remains unclear, and the current observations do not allow 
to resolve this controversial question. Moreover, given that 
other mechanisms such as competitive accretion (see sec- 
tion 4.3.) can explain massive star formation in the OHC, 
there is no need to invoke coalescence to build up massive 
stars. 
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